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Abstract 

Results for neutrino flavor oscillations and neutrino mixing mech- 
anisms, obtained from the analysis of the Sudbury Neutrino Obser- 
vatory (SNO), the SuperKamiokande (SK), CHOOZ, KamLAND and 
WMAP data, are used to calculate the effective neutrino mass < m„ > 
relevant for the neutrinoless double-beta decay (0v(3/3). The observ- 
ability of 0v(3(3 in the decay of 76 Ge is discussed within different light- 
neutrino mass hierarchies and by presenting a systematics on the avail- 
able nuclear matrix elements. 
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1 Introduction 



The knowledge about the properties of the neutrino has been dramatically 
advanced by various large-scale experiments, as recently reported by the 
SNO [1], SK [2], KamLAND [3], CHOOZ [4], and WMAP [5] collaborations. 
These experimental evidences have confirmed the existence of neutrino fla- 
vor oscillations and have set stringent limits to the neutrino mass-mixing 
mechanisms. A general overview of the latest experimental results is given in 
recent review articles by Valle [6], Bahcall et al. [7] and Elliot and Vogel [8]. 
Detailed discussions about the extracted values of the mixing angles, mixing 
amplitudes, and mass differences can be found in Refs.[9, 10, 11]. The im- 
plications of the latest results on the physics of electroweak interactions and 
dark-matter studies have been presented in Refs.[12]-[14]. 1 

In addition to the findings on neutrino flavor oscillations and the con- 
firmation of some of the theoretically predicted possibilities for the mixing 
and enhancement of the oscillations in the presence of matter [15], double- 
beta-decay experiments can provide complementary information on the na- 
ture of the neutrino and about its absolute mass scale [16]-[20]. This is a 
unique feature of the double-beta decay, which must be consistent with other 
scale- fixing measurements, like the WMAP measurements [5]. In the case of 
double-beta-decay measurements the knowledge about relevant nuclear ma- 
trix elements is crucial, as it is crucial to know the correct neutrino-mass 
spectrum for the analysis of the other type of measurements. The impli- 
cations of the results of the solar, atmospheric, reactor and astrophysical 
neutrino experiments upon double-beta-decay experiments have been stated 
already in several publications, see for instance [21]-[26]. To the wealth of pa- 
rameters involved in the analysis, like CP-phases, mixing angles and masses, 
one should add the nuclear-structure degree of freedom needed to extract the 
effective electron- neutrino mass [17]- [20]. 

At first glance, to physicists who are less familiar with nuclear-structure 
analysis, it may seem an easy task to produce the needed nuclear-structure 
information. Unfortunately it is not so because of several reasons: 

a) Double-beta-decay transitions take place in medium- and heavy-mass 
systems, where explicit shell-model calculations are unfeasible, unless severely 

1 Because of the large amount of publications in the field we focus our attention on the 
most recent ones, since most of the valuable previous literature has been quoted in the 
papers which we have included in the present list of references. 
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truncated valence spaces are used. 

b) The sensitive part of the calculations depends on the information about 
the structure of the states of the double-odd-mass nuclei. These intermediate 
states play an essential role in the second-order transition matrix elements 
entering the expression of the decay rate, and less is known about them, as 
compared with the relatively large amount of information gathered about 
double-even- mass nuclei. 

c) In dealing with medium- and heavy-mass nuclei one has to introduce 
approximations to obtain the participant wave functions and these approxi- 
mations are not unique, they vary from model to model. 

d) To assign a certain degree of significance to the already existing theo- 
retical results one has to define, first, what should be taken as the equivalent 
of the experimental confidence level, e.g. which models may be taken as ref- 
erences and what would be the confidence level assigned to them depending 
upon the used approximations. 

In the past, all of these features have been referred to as the uncertainties 
in the nuclear matrix elements and roughly estimated to be within factors 
of 2 to 3, with respect to some reference values. This aspect of the problem 
certainly deserves some attention, as we are going to discuss later on in this 
work, since there turns out to be a gap between the range of masses extracted 
from double-beta-decay studies, 0.4 eV to 1.7 eV [22], and those extracted 
from the other neutrino-related studies which yield upper limits of the order 
of 0.10 to 0.20 eV [24] or even lower [12]. There is a clear discrepancy between 
both sets of results concerning the observability of neutrinoless double-beta 
(Ov/3/3) decay. This issue has become a hot one, due to the recent claim [27] 
about the positive identification of neutrinoless-double-beta decay signals in 
the decay of 76 Ge (see also [28, 29]), from which a central value of the mass of 
the order of 0.39 eV was extracted [27]. We think that these aspects must be 
considered from points of view of both the neutrino physics and the nuclear- 
structure physics. To achieve this goal, in this work we discuss the constraints 
set by the oscillation and mass parameters on the effective neutrino mass 
measured in 0z//3/3-decay transitions and compare these neutrino masses with 
the ones coming from the analysis by using nuclear-structure information. 
We start from the best-fit mass-mixing matrix presented in [30] and from 
other estimates of the mixing matrix, i.e. the parametrization in terms of 
the mixing angle of solar neutrinos, and the estimation based on a maximum- 
mixing scheme [23]. 
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In the first part of the paper we review the basic elements of the theory 
and discuss the structure of the adopted neutrino mass-mixing matrix. We 
discuss a way to extract light-neutrino masses (m») from the observed mass 
differences and by combining them with the adopted neutrino mass-mixing 
matrix we calculate the effective neutrino mass relevant for the Ov/3/3 decay. 
In the second part of the paper we review the current nuclear structure 
information about the Ovflfl decay, by presenting the up-to-date values of the 
effective neutrino mass extracted from the adopted limits on the half-lives. 
In doing so, we have considered the range of variation for the nuclear matrix 
elements, calculated within definite classes of models. We have focused our 
attention to the case of the Ou/3/3 decay of 76 Ge. The analysis covers the 
calculated values of the nuclear matrix elements during the last fifteen years. 
This information is needed to estimate the plausibility of future double-beta- 
decay experiments. Finally, we discuss the observability of the Qvf3(3 decay 
in the context of the present results. 

2 Formalism 

2.1 Neutrino data 

Two- and three-generation analyses of neutrino data, provided by the so- 
lar and atmospheric observations and by the range of mass differences ex- 
plored in reactor-based experiments, have been performed by several groups 
[6]- [11]. The picture which emerges from these very detailed analyses of 
neutrino-flavor oscillations favours the Large Mixing Angle (LMA) solution 
of the Mikheyev-Smirnov-Wolsfenstein (MSW) mechanism [15]. Recently, 
the KamLAND collaboration [3] has confirmed the LMA solution and a cru- 
cial step towards the elucidation of the neutrino-mass hierarchies was given 
by the results of WMAP [5, 24], which fix a stringent upper limit for the 
scale of neutrino masses. A brief compilation of the adopted results is given 
in Table 1. As shown in this table, the SNO data are consistent with a value 
of the mass difference AmJ 2 of the order of 10~ 5 eV (solar-neutrino data), 
and another independent scale Am^ pa Aro^, of the order of 10~ 3 eV, has 
been determined from the analysis of the atmospheric-neutrino data, which 
is in the range of the sensitivity of the reactor-based measurements. Because 
of the independence of the determined mass differences, the global picture 
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is consistent with the existence of three active neutrino flavours. To these 
data, the information obtained by WMAP is adding the value of the upper 
limit of the sum of the three mass eigenvalues (light-neutrino masses only), 
which is of the order of 0.71 eV [5]. 

To calculate effective neutrino properties, like the effective electron-neutrino 
mass, < m v >, one needs to know the neutrino- mixing matrix U and the 
light-neutrino mass hierarchy (mi,m2,m 3 ) [16]. The determination of the 
matrix elements of U and the absolute values of the masses is the ultimate 
goal of any of the models of the neutrino and it is, of course, a matter of inten- 
sive effort. Out of the very rich, recently published list of articles dealing with 
the analysis of the SNO results, we have selected the results presented in the 
paper of Bandyopadhyay, Choubey, Goswami and Kar (BCGK) [30], together 
with the expression of the mixing matrix in terms of the solar-neutrino data, 
and the zeroth-order approximation of the mixing matrix assuming maxi- 
mum mixing, to perfom our calculations. Our choice is motivated by the 
fact that in the BCGK paper the best-fit value of U, with respect to the so- 
lar, atmospheric, and CHOOZ data, is written explicitly and the confidence 
level of the results is well established. 

The three-generation mixing matrix U can be written as 



This expression does not include CP violation. By performing a three- 
generation x 2 -analysis of the solar-neutrino and CHOOZ data, and by con- 
sidering the mass differences Am^ = Am^, Am^ ~ Am| = Am^ tm , the 
BCGK found that the best fit occurs in the LMA region with tan 2 d\^ ~ 0. 
This finding greatly simplifies the form of the mixing matrix U, because it 
narrows the value of U e s down to a very small range around U e s ~ [23, 30]. 
The best-fit form of U, reported in the BCGK paper, is 



The expression of the matrix U, considering U e 3 = and exploiting the 
solar and atmospheric mixing-angles data, reduces to (see Eq. (1)) 



/ C13C12 S12C13 S13 N 

U = — S12C23 — S23 s 13 c 12 c 23 c 12 — s 23 s 13 s 12 S23 c 13 
\ S23 s 12 — •S13C23C12 — S23 c 12 ~ Sl3 s 12C23 C23CI-3 / 



(1) 
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( C12 S12 ^ 

U = -S12C23 C23C12 S23 

V S23S12 -S23C12 c 23 / 

Finally, the matrix U, written in terms of maximum mixing (sin#i 2 
cos^23 = ^75), is of the form 



(3) 



U = 



1$ 

_ 1 

2 

v \ 



-4 \ 

1 

73/ 



1 

1 

2 
_ 1 

2 



(4) 



Only the first row is relevant for the electron-neutrino mass. The next 
step consists of the definition of a neutrino mass hierarchy, that is the relative 
order of the mass eigenvalues, which cannot be fixed only by the measured 
squared mass differences. 

In order to estimate the possible range of the we define the relative 
scales 

mi = fm 2 , m 2 = gm 3 (5) 
for the normal mass hierarchy (mi w m 2 < 777.3), and 



mi = fm 2 , m 3 = gm x 



(6) 



for the inverse [m\ m m 2 > m.3) and degenerate (777.3 ~ m 2 ~ mi) mass 
hierarchies. To these factors we have added the information related to the 
scale of the mass eigenvalues, which is determined by the extreme value 



mo = T 



(7) 



where the value of Q u is taken from the WMAP data (see Table 1). The fac- 
tors / and g are determined in such a way that the resulting masses mi(f, g) 
obey the observed mass differences, hereafter denoted as Am 2 (Am 31 ~ 
Am| 2 ) and Sm 2 (Am\ 2 )- We are restricted to light-neutrino masses, as said 
before. The numerical analysis was performed by assuming the above given 
scalings and by finding the values of (/, g) which are solutions of the equations 
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for the normal hierarchy case and 
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for the inverse and degenerate mass hierarchies. The use of the scale m fixes 
the limiting values of / and g at 
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(11) 



0<g < . = = (12) 



for the nearly degenerate masses. 

In the above expressions the factor r is given by the ratio between the 
solar and atmospheric squared mass differences 

Therefore, the variation of the parameters / and g is effectively restricted by 
the actual value of r and m . 
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For each set of allowed values of (/, g) and for each of the hierarchies 
considered we have calculated m*. The effective neutrino mass < m v >, 
relevant for the 0u(3(3 decay, is given by [17]- [20] 

3 

< m v >± = J^rriiXil U ei | 2 = m x Ul x ± m 2 U^ 2 , (14) 
i=i 

since for the adopted best fit U e3 [30]. We have consistently neglected 
CP violating phases and written A, = ±1, since the fit of [30] was performed 
under the assumption of CP conservation. In Table 2 we give, for each of the 
adopted forms of the mixing matrix U, the range of values of the calculated 
effective electron-neutrino masses. These values correspond to the limiting 
values of / and g, given in the previous equations (10)-(12). As can be seen 
from this table, the largest value which one can obtain for < m u > is of the 
order of 0.24 eV, and the smallest one is of the order of 0.7 xl0~ 4 eV, both 
for the degenerate mass hierarchy. Notice that the larger value is of the order 
of the mass scale extracted from the results of WMAP and it will certainly 
depend upon new results for VL V . A value of Q v < 0.5 eV [24] would then 
give a mass of the order of 0.16 eV, while the estimate Q v < 0.18 eV [12] will 
reduce it to the more stringent value of 0.06 eV. 

2.2 Nuclear matrix elements 

The implication of these results for < m v > upon the rates of 0v(3/3 decay is 
easily seen if one writes 2 the corresponding half-life, ty U 2 \ as 

(*')"' = {^)'°SS • < 15 > 

where the factor is defined as 

Ct2 = Gr\M§»\l- XF )) 2 , (16) 

in terms of the nuclear matrix elements, Mq^(1 — Xf), an d the phase-space 
factors, Gf v \ entering the mass term of the transition probability [17]. 

2 Only the mass sector of the half- life will be considered in the following analysis. Com- 
plete expression of the half-life, including right-handed currents, can be found in [17]. 
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There are several aspects concerning Eq. (15) which are worth of men- 
tioning: 

a) In the event of a successful measurement of Qv(3(3 decay and consider- 
ing the improving information emerging from neutrino-related measurements, 
Eq. (15) may be viewed as a crucial test for nuclear models, since the calcu- 
lated matrix elements reside in the factor C^^. 

b) If one assigns a certain confidence level to nuclear-structure calcula- 
tions, by fixing the value of C^, and takes the range of values of the effective 
neutrino mass extracted from neutrino-related measurements, Eq. (15) may 
be viewed as a criterium for determining the observability of 0uf3(3 decay. 

c) In the event of a positive measurement of Ov/3/3 decay and considering 
a reliable estimate of the nuclear matrix elements, Eq.(15) may be viewed as 
a consistency equation for the effective neutrino mass seen in double-beta- 
decay as compared with the one extracted from neutrino-related experiments. 

Let us start with the discussion of the nuclear-structure related infor- 
mation, contained in Cj^. The ultimate goal of nuclear-structure models 
is, in fact, the prediction of observables based on the knowledge about nu- 
clear wave functions down to the needed level of accuracy. In the case of 
Qv(3(3 decay studies, to achive this ultimate goal one needs to surpass sev- 
eral requirements, some of which are purely technical and some of which 
are conceptual. Among the technical barriers one has, of course, the unfea- 
sibility of large-scale shell-model calculations, prohibited by hardware con- 
straints. Among the conceptual requirements one has the realization that a 
prediction of a neutrinoless double-beta-decay rate should always be accom- 
panied by other model predictions, like single-beta-decay, electromagnetic 
and particle-transfer transitions involving the nuclei which participate in the 
double-beta-decay transition under consideration. We stress the point that, 
in our experience, the study should be conducted on the basis of a case-by- 
case analysis. 

Most of the current nuclear-structure approximations are based on the 
proton-neutron quasiparticle random phase approximation (pnQRPA). This 
is a framework where proton-neutron correlations are treated as basic build- 
ing blocks to describe the nuclear states which participate in a double-beta- 
decay transition. The pnQRPA formalism is rather well known and it has 
been discussed in a large number of publications during the last forty years. 
For the sake of brevity we are not going to present it here again, rather we 
would like to refer the interested reader to [17] for details. In particular, 
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the sensitivity of the pnQRPA method to values of specific parameters of 
the interactions, like the sensitivity to the renormalization of the particle- 
particle (proton-neutron) coupling, has been a matter of intensive studies. 
Again, we would like to refer to [17] for details concerning this point as well 
as concerning the large number of extensions of the pnQRPA method, their 
successes and failures. Restricting ourselves to a very elementary theoretical 
background, we can say that the standard procedure, applied in the literature 
to calculate the Ozz/^/S-decay rate, involves three major components: 

a) The calculation of the spectrum of the intermediate double-odd-mass 
nucleus with (A, N±l, Z^fl) nucleons. The pnQRPA is an approximate diag- 
onalization in the one-particle-one-hole, lp-lh, (or two-quasiparticle) space 
and it includes the effects of 2p-2h ground-state correlations by means of 
the backward-going amplitudes. Since the calculations are based on a quasi- 
particle mean field one forces the breaking of certain symmetries, like the 
particle-number symmetry by the use of the BCS approximation, and the 
isospin symmetry, by the use of effective proton and neutron single-particle 
states. The final results of the pnQRPA calculations will certainly be af- 
fected by these symmetry-breaking effects induced by the way in which we 
handle the nuclear interactions. Some attempts to cure for these effects 
have been implemented by means of enlarging the representation space, in- 
cluding Pauli-principle-related blockings and by performing self-consistent 
approaches beyond the quasiparticle mean field. As said before, the list of 
various extensions of the standard pnQRPA is too long to be commented 
here in detail. A fairly complete list of references about the set of extensions 
of the pnQRPA is given in [17]- [19]. We will generally refer to these approx- 
imations as pnQRPA-related ones. In this paper we shall show the results 
based on this family of approximations 3 . In addition, we quote the results 
of the available shell-model calculations. 

b) The calculation of the leptonic phase-space factors, as dictated by 
the second-order perturbative treatment of the electroweak interaction. At 
the level of the minimal extension of the Standard Model (SM) lagrangian 
(mass sector only), these phase-space factors can be easily calculated, and the 
values of them should be rather universal causing no source of discrepancies 

3 We shall explicitly quote the sources from where the results have been taken in order 
to avoid here a repetition of the details of each formalism, since the aim of the present 
section is not to present a critical review of the theories but rather to show their results 
to give an idea about the spread in the values of the relevant nuclear matrix elements. 
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in the calculations, except for the adopted value of the axial- vector coupling 
qa- At the level of the two-nucleon mechanism this value is currently fixed 
& t 9a = 1-254 but for the medium- heavy and heavy nuclei an effective value 
of gA = 1-0 has also been used. In this work we adopt the conservative 
estimate of qa = 1.254. Expressions for the phase-space factors, for theories 
beyond the minimal extension of the SM lagrangian, i.e. for left-right and 
right-right couplings, have been listed exhaustively in the literature (see e.g. 
[17]-[18]) and their values are well defined, too. In going beyond the two- 
nucleon mechanism one has to consider, also, the momentum dependence 
of the operators, which will reflect upon the structure of the phase-space 
factors. This is also true for the case of calculations where one is including 
p-electron wave effects and/or forbidden decays. 

c) The calculation of the matrix elements of the relevant current opera- 
tors which act upon the nucleons. These operators are also well known and 
their multipole structures are derived from the expansion of the electroweak 
current [17]. In the present calculation we have considered the standard type 
of operators, without introducing further momentum dependences in them, 
as originating from the electroweak decay at the quark level [31]. 

A compilation of the values of nuclear matrix elements and phase-space 
factors can be found in [17]. The current information about the status of 
0uf3f3 decays is reported in [32]- [36]. 

Tables 3 and 4 show the set of double-beta-decay systems where experi- 
mental searches for signals of the Qv(3(3 are conducted at present or planned 
for the next generation of double-beta-decay experiments. The tables con- 
tain the experimental lower limits for the 0u(3f3 half-life [37]- [46], the full 
range of variation of the nuclear matrix elements, as contained in the factors 
and as they are predicted by different models [17], the values of the 
model-dependent factor [17] 

F N = t$C££ = (< m„ > /m e y 2 , (17) 

the calculated phase-space factors Gf u \ and the extracted values of the 
effective neutrino mass. 

In Table 4 only a sub-group of calculations are presented, namely the 
ones based on the plain spherical pnQRPA approach of [17] (third column). 
These results are compared with our present calculations shown in the fourth 
column. 
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In the following, some brief details about the present pnQRPA calcula- 
tions are given. They have been done by following the procedure outlined 
in [17]. The two-body nuclear interactions were constructed by using the G- 
matrix interaction of the Bonn type including two to three major harmonic- 
oscillator shells around the proton and neutron Fermi surfaces. The spherical 
Woods-Saxon potential was used to generate the single-particle energies and 
small adjustments of these energies were done in the vicinity of the Fermi 
surfaces to reproduce the low-energy quasiparticle spectra of the neighbour- 
ing odd-mass nuclei. Following the criteria which we have advanced above, 
the various parameters involved in the calculations have been fixed by repro- 
ducing the known data on single-beta-decay transitions around the nuclei of 
interest for the double-beta-decay transitions which we are considering here. 
No further adjustement of the proton- neutron particle-particle coupling con- 
stant [17] is introduced once the known single-beta decay observables are 
reasonably reproduced. 

As one can see, our present results are in good agreement with the other 
pnQRPA calculations, except for 136 Xe where our calculation gives a larger 
matrix element than the other calculations. This deviation might occur due 
to the semi-magic nature (the neutron shell is closed) of 136 Xe, forcing the 
transition from the two-quasiparticle description to the particle-hole descrip- 
tion. 

If one compares the extracted neutrino masses of Table 3 with the ones 
given in the previous section, it becomes evident that the present gener- 
ation of experiments is rather insensitive to the effective neutrino 
mass coming from the best fit of the solar+atmospheric+reactor data, ex- 
cept for the Heidelberg-Moscow experiment if one takes the range of values 
(< m v >=0.11 eV-0.56 eV) reported in [27]. If one takes the heaviest possible 
effective mass < m v >^ 0.24 eV, which is favoured by the inverse and degen- 
erate mass hierarchies (see Table 2), one sees that it is outside the range of 
the present limits fixed by double-beta-decay experiments, with the possible 
exception of the decay of 76 Ge, which just barely reaches this estimate. Nat- 
urally, the values of < m v > extracted from the experimental lower limits of 
tfj^ are model dependent, since the connection between the half-lives and the 
effective neutrino mass is given by the nuclear- model-dependent factors F N . 
As a reference value, for < m v >^ 0.24 eV one obtains = 4.53 x 10 12 (see 
Eq.(17)) which is to be compared with the estimate (see Table 3, case of 76 Ge 
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) F N = 3.55 x 10 11 -> 7.20 x 10 12 , computed by assuming 4/2 > 2 - 5 x 1q25 J t 
[29] and taking into account the total span of the calculated nuclear matrix 
elements. 

With reference to the results shown in Table 4, the span in the effective 
neutrino mass is smaller (see the last column of Table 4) than when all the 
available model calculations are included (see the last column of Table 3). 
For the case of 76 Ge the spherical pnQRPA gives a span of < m u >=0.30 eV 
-0.43 eV in the effective mass. This means that to reach the neutrino-mass 
value resulting from the neutrino data, one definitely needs larger matrix 
elements than the ones produced thus far by the spherical pnQRPA model, 
and/or longer half-lifes than the present measured limits. These observations 
will be discussed more in the next section. 



2.3 pnQRPA Matrix elements for 76 Ge 

Table 5 shows the results of the matrix elements, corresponding to the mass 
sector of the neutrinoless double-beta decay in 76 Ge, calculated within the 
family of the pnQRPA-related models [47]- [54]. The standard spherical pn- 
QRPA method gives results which are of the order of C^ m w 5 - 8 x 10~ 14 
in units of yr _1 , with the exception of the result presented in [54], which 
yields to a magnitude of the order of 1.85 x 10~ 14 yr _1 , and the one of [47] 
where the pnQRPA value is 1.12 x 10~ 13 yr _1 . These factors translate into 
the ranges of the nuclear matrix elements 4 and effective neutrino masses 
which were shown, previously, in Tables 3 and 4. The results of the other, 
pnQRPA-related, approximations seem to be less stable and they deviate 
more from the central range of C$ m ~ 5 — 8 x 10~ 14 yr _1 . In analyzing the 
results of [54] one can notice that the largest value does not differ much from 
the standard pnQRPA value, although is has been obtained by using a more 
involved approximation. By using the phase-space factors listed in Table 
4, we arrive at the central value for the matrix elements in the pnQRPA, 
namely 

Mgx ) (l -XF) pnQRPA = 3.65 . (18) 

4 Notice that the results of [47], which are relevant for the analysis performed in [27], 
are only 1.3 times larger than the average pnQRPA matrix element. 
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The corresponding value for the latest large-scale shell-model calculation [57] 
is given by 

Mg?(l-X F ) shcll _ modcl = 1.74 . (19) 

Therefore, the latest shell-model results [57] and the centroid of the pnQRPA 
results differ by a factor of the order of 2. In terms of the effective neutrino 
mass, using the half-life limit 2.5 x 10 25 yr, recommended in [29], these matrix 
elements lead to 

< rn„ > pnQRPA = 0.35 eV , (20) 
for the pnQRPA estimate, and 

< m v > sh cii-modci = 0.74 eV , (21) 

for the shell- model estimate of the matrix element. It means that to go to 
masses of the order of 0.24 eV, as required by WMAP, one needs larger nu- 
clear matrix elements than the ones given by the pnQRPA or by the available 
shell-model results. In fact, to reach the WMAP limit one would need the 
value 

MS ) (l-X F ) cxp en m e„tal= 5.36, (22) 

which is ~ times larger than the reference pnQRPA value given in (18). 
The largest matrix element listed in Table 5, coming from the VAMPIR 
approach [55], would yield to the value < m v >vampir — 0.19 eV, which 
just touches the value < m v > = 0.24 eV coming from the analysis of the 
neutrino-related data. However, it is appropriate to point out here that the 
VAMPIR matrix element is to be considered unrealistically large because in 
the calculations of [55] no proton-neutron residual interaction was included. 
Finally, our present value 

M^ ) (1-XfC C qrpa = 3.33 (23) 

(see Table 4) is consistent with the central value (18), and it yields an effective 
neutrino mass 

< m v >£3£ A = 0.39 eV (24) 
if one takes for the half-life the lower limit recommended in [29], and 

< m v >£5£ A = °- 50 eV ( 25 ) 

if one takes for the half-life the value 1.5 x 10 25 yr given by the Heidelberg- 
Moscow collaboration [27]. 
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3 Observability of the neutrinoless double beta 
decay 

To grasp an idea about the observability of the Ovflfl decay in other systems, 
we can compare the values of F N , of Table 3, with the ones obtained by using 
an effective neutrino mass of 0.24 eV, corresponding to F N = 4.53 x 10 12 . 

Figure 1 shows the comparison between the values of Fn of Eq.(18), listed 
in Table 3, and the values corresponding to the effective neutrino masses 
< m v >= 0.24 eV and 0.39 eV. The interval between upper and lower values, 
for each case, represents the span of the calculated nuclear matrix elements. 
For the case of 76 Ge the prominent upper value is given by the unrealisticaly 
large nuclear matrix element of [55]. 

The results shown in Figure 1 indicate a departure with respect to the ex- 
perimental limits by orders of magnitude, excluding the case of 76 Ge which is 
closer but still outside of the range consistent with the solar+atmospheric+reator 
neutrino data. 

Thus the issue about the observability of the Ovfifi decay relies upon the 
estimates for the effective neutrino mass and upon the estimates of the rel- 
evant nuclear matrix elements. While in some cases the differences between 
the calculated matrix elements are within factors of the order of 3, in some 
other cases the differences are much larger. It shows one of the essential fea- 
tures of the nuclear double-beta decay, namely that case-by-case theoretical 
studies are needed instead of a global one [17]. 

4 Conclusions 

To conclude, in this paper we have presented results on the effective neu- 
trino mass, as obtained from the best-fit mass-mixing matrix U determined 
from the analysis of solar+atmospheric+reactor+satellite data, and com- 
pared them with the values extracted from neutrinoless double-beta-decay 
experiments. 

The effective electron-neutrino mass extracted from the neutrino-related 
experiments, < m u >m 0.24 eV, does not compare with the central value 
of < m u >m 0.39 eV , reported in [27] and obtained by using the nuclear 
matrix elements calculated in [47]. It does not compare, either, with the 
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values given by the standard pnQRPA model, after taking into account the 
span in the calculated matrix elements. 

To explain for the difference between the above results, we have compiled 
a systematics of the calculated nuclear matrix elements and performed addi- 
tional pnQRPA calculations. In the case of 76 Ge, and if one adopts for the 
half-life the upper limit of 2.5 x 10 25 yr suggested in [29], the nuclear matrix 
elements needed to yield the desired effective neutrino mass are larger than 
any of the known nuclear matrix elements calculated in the framework of the 
spherical pnQRPA. This is conclusion also holds for the available shell-model 
results. 

The present knowledge of the involved nuclear matrix elements shows that 
the sensitivity of the Qv/3/3 experiments is still far from the estimate coming 
from neutrino oscillations. However, the needed sensitivity is potentially 
achievable by the next generation of experiments. 
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Figure Captions 

Figure 1: Factors F/v, of Eq. (18), for each of the systems of Table 3. 
The lines are drawn to guide the eye. The interval between the upper and 
lower lines, for each case, represents the span of the calculated nuclear matrix 
elements. The results corresponding to < m v >= 0.24 eV and < m v >= 
0.39 eV are shown as horizontal lines. 
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Table 1: Current limits on neutrino-mass differences. The values listed are 
a compilation of the results from the SNO [1], SK [2], KamLAND [3] and 
WMAP [5]. 

5m 2 12 = 8m 2 solar 5 x 10" 5 eV^ ^ 1.1 x 10" 4 eV 2 

5m 2 23 = 5m\ tm 10~ 3 eV 2 -> 5 x 10" 3 eV 2 

sin 2 2^ so i ar « 0.86 

sin 2 2# atm w 1.0 

ft„ < 0.71 eV 



Table 2: Calculated effective electron-neutrino masses < m v > ± . Indicated 
in the table are the mass hierarchy and the adopted mixing matrix. The 
values are given in units of eV. The results listed as extreme have been 
obtained by using the extreme upper values of / and g of Eqs.(10)-(12). 
The adopted values for the mass differences are 5m 2 i 2 = 7.1 x 10~ 5 eV 2 , 
<5m 2 23 = 2.7 x 10~ 3 eV 2 , and m =0.24 eV. The mixing matrix U(a) is taken 
from the best fit of [30], U(b) is based on the largest values of the solar and 
atmospheric mixing angles, and U(c) assumes maximum mixing. 



Hierarchy 




< m t 




U(a) 


U(b) 


U(c) 




Normal 


(mi = 0) 


< m v 


> 


-0.010 


-0.012 


-0.019 








< m v 


>+ 


0.011 


0.012 


0.019 






(extreme) 


< m v 


>_ 


0.105 


0.086 


-0.769 xHT 


-4 






< m v 


>+ 


0.231 


0.231 


0.231 




Inverse 


(m 3 = 0) 


< m v 


>_ 


0.105 


0.087 


-0.153 xlO" 


-2 






< m v 


>+ 


0.234 


0.235 


0.235 






(extreme) 


< m v 


>_ 


0.108 


0.088 


-0.749 xlO" 


-4 






< m v 


>+ 


0.237 


0.237 


0.237 




Degenerate 


(extreme) 


< m v 


>_ 


0.107 


0.088 


-0.715 xlO^ 


-4 






< m u 


>+ 


0.237 


0.237 


0.237 
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Table 3: 0u(3(3: Model-dependent estimates and experimental limits. The 
double-beta-decay systems are given in the first column. The factors 
are given in units of yr _1 and their values are shown within the intervals 
predicted by different nuclear-structure models, like the Shell Model (a), the 
Quasiparticle Random Phase Approximation (b), the pseudo SU(3) model 
(c), and various other models (d). The value = 1.254 is used. The 
quantities ty U J are the experimental lower limits of the half-lives, in units of 
years. The corresponding references are quoted in brackets. The factor Fo- 
llower limit) is shown in the fourth column and the values are given within 
the intervals provided by the factors C^fy . The last column shows the range 
of variation of the extracted effective neutrino mass (upper limits) in units 
of eV. The coefficients are taken from [17], except for the case of 124 Sn 



[58]. 


System 










AOu) 
l l/2 




Fn 


< m v > 


48 Ca 


(1.55-4.91) 10" 


-14 


(a) 


9.5 


10 21 


[37] 


(1.47-4.66) 10 s 


(23.7-42.1) 




(9.35-363) 10- 


-15 


(b) 








(8.88-345) 10 7 


(8.70-54.2) 


76 Ge 


(1.42-28.8) 10" 


-14 


(d) 


2.5 


10 25 


[29] 


(3.55-72.0) 10 11 


(0.19-0.86) 


82 Se 


(9.38-43.3) 10- 


-14 


(d) 


2.7 


10 22 


[39] 


(2.53-11.7) 10 9 


(4.73-10.2) 


96 Zr 


(9.48-428)10- 


15 


(b) 


1.0 


10 21 


[40] 


(9.48-428) 10 6 


(24.7-166) 


100 Mo 


(0.07-2490) 10 


-15 


(b) 


5.5 


10 22 


[41] 


(0.38-13700) 10 7 


(1.38-262) 


116 Cd 


(5.57-66.1) 10- 


-14 


(b) 


7.0 


10 22 


[42] 


(3.90-46.3) 10 9 


(2.37-8.18) 
(1.38-2.18) 10 3 


124 Sn 


(2.29-5.70) 10- 


-13 


(b) 


2.4 


10 17 


[43] 


(5.50-13.7) 10 4 


128 Te 


(1.71-33.6) 10" 


-15 


(b) 


8.6 


10 22 


[44] 


(1.47-28.9) 10 8 


(9.51-42.1) 


130 Te 


(1.24-5.34) 10- 


-13 


(b) 


1.4 


1Q 23 


[44] 


(1.74-7.48) 10 10 


(1.87-3.87) 


136 Xe 


(2.48-15.7) 10- 


14 , 


(a,b) 


4.4 


1Q 23 


[45] 


(1.09-6.91) 10 10 


(1.94-4.89) 


150 Nd 


(4.78-77.4) 10" 


13 


(b,c) 


1.7 


10 21 


[46] 


(8.13-132) 10 8 


(4.45-17.9) 
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Table 4: Calculated phase-space factors G\ u> and calculated nuclear matrix 
elements, using the formalism of the spherical pnQRPA, for some of the 
double-beta emitters included in Table 3. The phase space factors are given 
in units of yr -1 and the dimensionless matrix elements are scaled by the 
nuclear radius [17]. The third column, indicated as N.M.E., gives the extreme 
values of the nuclear matrix elements Mq' t (1 — xf) reported in the literature 
(see the captions to Table 3), and the fourth column, indicated as N.M.E. 
(this work), gives the results of the present calculations for Mq T (1 — Xf)- 
The last column shows the range of values for the effective neutrino mass, in 
units of eV, extracted from the results given in the third and fourth columns. 



System 


Gf v) x 10 14 


N.M.E. 


N.M.E. (this work) 


< m v > 


4 «Ca 


6.43 


1.08-2.38 




8.70-19.0 


76 Ge 


0.63 


2.98-4.33 


3.33 


0.30-0.43 


82 Se 


2.73 


2.53-3.98 


3.44 


4.73-7.44 


96 Zr 


5.70 


2.74 


3.55 


19.1-24.7 


100 Mo 


11.30 


0.77-4.67 


2.97 


1.38-8.42 


116 Cd 


4.68 


1.09-3.46 


3.75 


2.37-8.18 


128 Te 


0.16 


2.51-4.58 




9.51-17.4 


130 Te 


4.14 


2.10-3.59 


3.49 


1.87-3.20 


136 Xe 


4.37 


1.61-1.90 


4.64 


0.79-2.29 
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Table 5: Calculated nuclear matrix elements for the case of Ge. The val- 
ues C$ m are given in units of yr -1 . The adopted value for the half-life is the 
value recommended in [29], = 2.5 x 10 25 yr. Indicated in the table are the 
models used to calculate the nuclear matrix elements, which are taken from 
the references quoted in the last row of the table. The abreviations stand for 
the proton-neutron quasiparticle random-phase approximation (pnQRPA), 
particle-number-projected pnQRPA (pnQRPA (proj.)), proton-neutron pair- 
ing pnQRPA (pnQRPA+pn pairing), the renormalized pnQRPA (RQRPA), 
the second pnQRPA (SQRPA), the self-consistent renormalized pnQRPA 
(SCRQRPA), the fully renormalized pnQRPA (full-RQRPA), and the vari- 
ation after projection mean-field approach (VAMPIR). The model assump- 
tions underlying these theories are presented in the quoted references. 



C(0) 

mm 


F N x 10~ 12 


Theory 


Reference 


1.12 x 10- 13 


2.80 


pnQRPA 


[47, 48] 


6.97 x 1(T 14 


1.74 


pnQRPA 


[31] 


7.51 x 1(T 14 


1.88 


pnQRPA (proj.) 


[31] 


7.33 x 1(T 14 


1.83 


pnQRPA 


[49] 


1.42 x 1(T 14 


0.35 


pnQRPA+ pn pairing 


[49] 


1.18 x 1CT 13 


2.95 


pnQRPA 


[50] 


8.27 x 1(T 14 


2.07 


pnQRPA 


[51] 


2.11 x 1(T 13 


5.27 


RQRPA 


[52] 


6.19 x 10" 14 


1.55 


RQRPA+ q-dep. operators 


[52] 


1.8-2.2 x 1(T 14 


0.45-0.55 


pnQRPA 


[53] 


5.5-6.3 x 1(T 14 


1.37-1.57 


RQRPA 


[53] 


2.7-3.2 x 1(T 15 


0.07-0.08 


SCRQRPA 


[53] 


1.85 x 1(T 14 


0.46 


pnQRPA 


[54] 


1.21 x 1(T 14 


0.30 


RQRPA 


[54] 


3.63 x 1(T 14 


0.91 


full-RQRPA 


[54] 


6.50 x 1(T 14 


1.62 


SQRPA 


[54] 


2.88 x 1CT 13 


7.20 


VAMPIR 


[55] 


1.58 x 1(T 13 


3.95 


Shell Model 


[56] 


1.90 x 10~ 14 


0.47 


Shell Model 


[57] 
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